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ACTIVATION, GROWTH, AND DETACHMENT OF BOILING BUBBLES IN WATER FROM 
ARTIFICIAL NUCLEATION SITES OF KNOWN GEOMETRY AND SIZE 
by John R. Howell and Robert Siegel 
Lewis Research Center 

SUMMARY 

Artificial nucleation sites in the range from 0.0015 to 0. 032 inch in radius were 
drilled in a polished metal surface by either electron-beam or mechanical drilling. The 
formation of boiling bubbles was studied by electrically heating the surface in saturated 
distilled water and by obtaining surface temperature, heat flux, and photographic data. 
The surface -temperature elevation above saturation required to produce ebullition from a 
pre-existing vapor nucleus is given as a function of site radius. The results are com- 
pared with several available analyses including one proposed in the present report. The 
analyses do not compare well with the data, thus a basic feature is still missing from 
the analytical models. In general, the analyses predict that bubbles will start to grow at 
surface temperatures lower than those measured experimentally. Data giving bubble 
growth as a function of site radius and surface temperature elevation are provided for 
saturated water and are compared with theoretical growth relations. At bubble detach- 
ment in saturated liquid, the buoyancy and surface-tension forces are found to be 
approximately equal. 


INTRODUCTION 

The boiling phenomena have been the subject of numerous research efforts as 
surveyed, for example, in references 1 and 2. These investigations have revealed that 
many factors affect boiling, and, at present, there is no correlation for heat -transfer 
coefficient that satisfactorily accounts for all parameters involved. As a consequence, 
some recent research has been undertaken with the view that a quantitative understanding 
of the boiling problem will be achieved only by a return to very basic studies of the 
boiling process, starting with the formation of individual bubbles and their growth and 
departure. The conditions under which bubbles will be produced from a specified cavity 
are the principal subject of the present report. 



Analytical predictions of the criteria for the initiation of bubble growth have been 
made in references 3 to 7, but only limited experimental data are available (refs. 5 
and 8 to 11) that relate bubble formation to a known site size and geometry. A satis- 
factory verification of these analyses has therefore not been possible. The analyses of 
bubble incipience criteria by Hsu (ref. 3) and by Han and Griffith (ref. 4) are quite 
similar. Both are based on the assumption that a surface cavity has a hemispherical 
vapor cap over it that serves as a bubble nucleus, and both assume that this nucleus 
will begin to grow only when the thermal layer of superheated liquid adjacent to the sur- 
face is of a sufficient thickness. The chief difference in the two analyses is in the 
assumption of how large the thermal-layer thickness must be to produce growth for a 
vapor nucleus of a given height. With an assumed relation between the thermal-layer 
thickness and the bubble nucleus size at incipience, a relation between active cavity size 
and surface -temperature elevation is derived. The transient conduction equation is 
used to provide the temperature distribution in the thermal layer, and the required tem- 
perature difference for bubble growth is found from the Clausius -Clapeyron equation. 

In reference 3, the range of active cavity size is given for a known surface -temperature 
elevation above saturation or for a specified constant heat flux at the wall. Reference 4 
discusses the case of specified surface temperature. 

Griffith and Wallis (ref. 5), among others, used the Clausius -Clapeyron equation 
and the Gaussian expression relating surface tension and radii of curvature to the pres- 
sure difference across the bubble interface to determine the active cavity radius. In 
contrast with the results of references 3 and 4, this relation does not indicate a maxi- 
mum cavity size and predicts that only small values of the temperature difference 
between the wall and the saturated liquid are required to initiate bubble growth as the 
cavity radius becomes large. The derivation is based on having a liquid of uniform 
temperature surrounding the bubble, and the liquid must be superheated for the bubble 
to grow. The relations given in references 3 and 4 do not have the restriction of a 
uniform liquid temperature. 

For a heated surface in a saturated bulk liquid, Madejski (ref. 6) includes the 
effect of the temperature gradient around the bubble causing a variation of bubble curva- 
ture. The wall- to saturated-liquid-temperature difference at incipience is found to be 
that predicted in reference 5 multiplied by a factor larger than one. For the conditions 
of the present experiments, the factor produced negligible change from the results of 
reference 5. 

Leont'ev and Kirdyashkin (ref. 7) provide an analysis based on free -convection 
theory for predicting the thickness of a thermal layer on a horizontal surface. The 
thickness is predicted to depend on the temperature difference between surface and bulk 
fluid to the minus one-third power. This thermal-layer thickness is then inserted into 
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the incipience analysis of Hsu (ref. 3) to predict the dependence of active cavity radius 
upon temperature. 

The actual trend with cavity radius of wall- to saturated-liquid-temperature differ- 
ence at incipience has not been completely resolved experimentally. As shown in refer- 
ences 5 and 11, the predicted incipience temperature differences are generally lower 
than those observed experimentally. In the work reported in reference 5, cavities 
0. 0027 inch in diameter were pricked in a surface with a sharpened needle. The cavities 
produced bubbles at a surface superheat of 20° F compared with a value of 3° F as 
predicted by the analysis of reference 5. In reference 11, cavities with diameters from 
0. 001 to 0. 005 inch were formed by etching with acid. For water at atmospheric pres- 
sure, the measured surface superheats required to provide bubbles were about 14° F 
compared with their predicted values of 2° to 4° F. 

Since the activation theories postulate that a vapor nucleus is initially present, they 
apply only to a range of cavity sizes in which vapor (or gas) is entrapped. The maximum 
cavity diameter thus should be limited to the value which will not allow replacement of 
vapor within the cavity by the surrounding liquid. Bankoff (ref. 12) presents a method by 
which this size may be established for certain systems. 

The present study was undertaken to determine experimentally the surface- 
temperature elevation above saturation required for a cavity to become active. Because 
of the many parameters which could affect the system, a simple configuration was 
chosen for the heated surface and cavity. Yatabe and Westwater (ref. 10) have indicated 
that cavity shape appears to have little effect on nucleating characteristics; thus, the 
only cavity parameter varied systematically in the present study was the diameter. The 
only system parameter varied was surface temperature (and therefore, of course, 
surface heat flux). The water was at saturation for all the tests. In addition to cavity 
activation data, high-speed motion pictures were taken at and near incipience of bubble 
formation to study the nature of the vapor nucleus and the bubble growth and departure 
characteristics from sites of known size and geometry. 

An earlier version of this report is given by reference 13. 


APPARATUS 

Heated Test Sections 

The test specimens were made from number 410 stainless -steel strips, 4 inches 
long, 1/2 inch wide, and 1/16 inch thick. This alloy was chosen for its magnetic 
properties which allowed fine polishing on available equipment that utilized magnetic 
clamps. The strips were polished on all edges and on one face to achieve a finish with 
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root mean square roughness less than 4 microinches. After polishing, a 20-gage Chromel- 
Alumel thermocouple was spotwelded to the center of the unfinished face of the specimen. 
The thermocouple leads were run parallel to the strip to reduce conduction losses. The 
strip was then placed in a mold which left only the finished face exposed, and epoxy 
resin was cast around the remaining sides providing a block of thermal insulation 
4 by 1 by 1/2 inches. The artificial nucleation sites were drilled into the polished face, 
and the entire specimen was cleaned ultrasonically in an alcohol bath and finally in a 
distilled -water bath. 


Artificial Sites 

Two types of artificial sites were studied. The first type was drilled by an 
electron-beam disintegration technique in the hope of producing extremely small diame- 
ters. However, the smallest site diameter attained by this method was 0. 0025 inch. 

In addition, the disintegration technique tended to leave a cratered area outside the 
actual hole diameter. This technique left some doubt as to the exact geometry of the 
cavity. The cavity depths ranged from 0. 003 to 0. 019 inch. A photograph of a typical 
cavity is shown in figure 1(a). 

The second type of artificial site was produced by a mechanical drill for the smaller 
holes and by an end mill for the larger diameters. These methods produced holes that 
were quite round at the surface and with no cratering evident (figs. 1(b) and (c)). Diam- 
eters near a minimum of 0. 003 inch were attained by this method using drills of 0. 002- 
inch nominal diameter. The largest diameter tested was 0. 064 inch. The cavity depths 
were from 0. 010 to 0. 020 inch. 

In most of the test strips, two sites of the same nominal diameter were drilled 
about 1 inch apart. Three strips were tested where either four or six holes of various 
diameters were drilled about 3/4 inch apart in a staggered pattern into the same strip; 
thus, the effect of site diameter could be observed while all the sites were under identi- 
cal conditions. 

The diameters of the cavities were measured using a metallograph, an optical 
device with a calibrated eyepiece, and the accuracy was estimated as ±0. 0002 inch for 
mechanically drilled sites. The photographs in figure 1, obtained using the metallo- 
graph, give a representation of the view seen through the calibrated eyepiece. Accuracy 
of the measurements for the electron-beam sites was about ±0. 0005 inch; the error was 
greater than for the mechanically drilled sites because of the uncertainty introduced by 
the cratering. 

Cavity depths were measured during the drilling process. A depth gage on the 
precision drill press provided this measurement for the mechanically drilled cavities 
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while the expended energy and time of disintegration were used to estimate the depths of 
the holes obtained by electron -beam disintegration. 


Pool Boiler 

The boiler consisted of a cylindrical stainless -steel tank, 14 inches in diameter and 
30 inches high (fig. 2). Two 8 -inch -diameter quartz windows were provided diamet- 
rically opposite each other for viewing or photography. Four 1000-watt heaters near 
the bottom of the tank provided auxiliary heat to maintain saturation conditions within the 
boiler. Two Chromel-Alumel thermocouples were immersed in the boiling fluid to 
monitor the bulk temperature. Power was provided through copper electrodes extending 
down through the boiling water. The test specimen was clamped to the electrodes with 
the test surface horizontal. The outer surface of the boiler, including the bottom, was 
wrapped with 2 inches of felt insulation. 


Auxiliary Equipment 

Power was supplied from a motor-controlled variable output alternating-current 
transformer. The alternating current provided to the test specimen produced no 
problems in the temperature measuring system. Calculations based on the analysis of 
reference 14 showed that the strip geometry and physical properties were such that 
temperature fluctuations resulting from the 60-cycle current were negligible. 

A high-speed camera with the associated timing equipment and lighting was available 
for photographing individual sites. Lenses and extension tubes were so arranged that 
the individual bubbles were magnified to practically fill the frame of the film. 


Instrumentation 

The current and the voltage across the test section were read with 1/4 percent 
alternating-current meters. The bulk water and test specimen thermocouple voltages 
were read on a laboratory -type potentiometer. The thermocouples were referenced 
to a 150°±l/4° F commercial thermocouple reference oven. 


PROCEDURE 

The boiler tank was drained and cleaned after every second day of experimentation. 
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An ultrasonically cleaned test specimen was clamped to the electrodes, and the tank 
was raised into position around the specimen. The boiler was then filled with distilled 
water, and the auxiliary heaters were turned on. The water was heated to saturation 
and was allowed to boil for a few hours to drive off any absorbed gases. The power to 
the auxiliary heaters was decreased to a level sufficient to maintain saturation while 
reducing the convective eddies formed by the boiling from these heaters. 

Power at a low level was then applied to the specimen, and the bulk temperature, 
current and voltage across the specimen, and specimen thermocouple readings were 
taken. Note was made as to whether the artificial site or sites were active, and the 
level of activity was observed (none, sporadic, intermittent, steady, etc.). The power 
was then increased several percent, and new readings were obtained. This process was 
continued until the artificial sites were producing bubbles continuously (except for a few 
cases where no nucleation occurred from sites that were previously active, possibly 
because the sites became filled with water). The power level was then reduced in steps 
of several percent each until all activity ceased. 

Photographs were taken when desired as the tests were too numerous to photograph 
every test condition. 

The surface temperature of the specimen was determined when the activity of the 
site was steady, was intermittent, or ceased, while the power was either increased or 
decreased. Thermocouple readings, which were taken on the back (insulated) side of 
the strip, were used to determine the temperature at the surface by assuming a one- 
dimensional heat -conduction process through the specimen thickness. For a plate of 
thickness a, insulated on one side and having uniform internal heat generation, the 
temperature difference between the two sides is given by 


"w 


T = SS. 
°“2k 


( 1 ) 


where T^ for the present experiments corresponds to the value provided by the speci- 
men thermocouple. (All symbols are defined in appendix A. ) 


RESULTS 

Activation of Nucleation Sites and Incipience of Boiling 

Definition of incipience point . - As shown by Hsu in the discussion appended to 
reference 13, a distinction must be made between the conditions necessary to cause a 
bubble nucleus to grow and the more stringent conditions necessary not only for growth 
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of the nucleus but also for detachment of the bubble from the site. The former will be 
termed "site activation" and the latter "incipience of boiling. " 

For a saturated system, the difference between the activation and incipience condi- 
tions will be small, because once a nucleus begins to grow it tends to continue its growth 
(if the thermal layer thickness does not decrease) until detachment occurs. For a sub- 
cooled condition, however, incipience may require a significantly higher surface tem- 
perature elevation above Tg^^ than that required for bubble activation. Activation 
thus occurs at the lower limit of conditions necessary for incipience. 

The word incipience has been defined in a number of other ways when applied to 
boiling. One common definition is that incipience corresponds to the surface- to bulk- 
temperature difference at which the curve of AT plotted against the heat flux for the 
entire surface begins to deviate from the free convection curve. Such a definition has 
little value in studying bubble production from individual sites involving isolated locations 
on the surface. 


Incipience can also be defined as the condition when the waiting time between depar- 
ture and initiation of a new bubble is not infinitely long. For the present experiments, 
however, this definition also has little meaning as each departing bubble left a vapor 
nucleus behind, and hence the waiting time was always zero. 

Hsu (ref. 3) and Han and Griffith (ref. 4) in their analyses consider site activation; 
that is, the conditions and temperature difference at which a given bubble nucleus begins 
to grow. This criterion is more valid for the lower limit of incipience in the experi- 
ments reported herein, but it still has some inherent difficulties. As the test- section 
surface temperature was increased, a point was reached where a bubble would grow at 
the nucleation site and then remain attached to the surface. This condition is one where 


the vapor nucleus grew and hence was active, but it still did not appear to be incipient 
boiling as a succession of bubbles was not produced from the site. At a slightly higher 
temperature level, the bubble detached, and a new bubble began to grow; but this bubble 
might linger at the surface for many seconds before detaching. For a larger test-surface 
temperature, the production of bubbles became more regular. Finally, at a still larger 
temperature difference, bubble production was so regular that a steady stream of bubbles 
appeared from the cavity. This leads to the following question: Which of these bubble- 
formation conditions and its corresponding surface- to saturation -temperature difference 
corresponds to incipience ? 

Because the irregularities in the nucleation phenomena may be due to the statistical 
fluctuations inherent in the thermal boundary-layer thickness, it could be conjectured that 
steady or regular bubble production could only occur when the superheated thermal-layer 
thickness reached a mean value which was greater than that used in the activation 
criteria of references 3 and 4. 


In the present study, the incipience point was examined with regard to two definitiona 
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One definition that provided an upper limit was that incipience corresponds to the value 
of AT at which regular bubble production occurred from a given site. It is obvious that 
such a definition involves possible errors in the judgement of the observer and may 
produce (AT)^ values above the theoretical site activation predictions. It also involves 
the experimental difficulty of defining when "regular” production is occurring from 
sites of differing size. If a large and a small site are compared in the region of in- 
cipience, the frequency of bubble emission is lower for the larger site as this site is 
producing larger bubbles; however, the percentage variations around a given mean 
bubble frequency seemed smaller for large sites. Thus, a small site producing bubbles 
might not be incipient according to this definition because of a lack of bubble regularity, 
while a large site producing bubbles at a lower frequency with bubble regularity would 
meet this incipience criterion. 

The second definition of incipience examined experimentally was that it corresponds 
to the minimum AT at which bubbles still grow and detach from the surface. Below 
this AT a bubble might still be present at a site, but it would remain attached to the 
surface. This criterion would be expected to correspond more closely with a theory 
predicting site activation but with no consideration as to whether a stream of bubbles is 
produced. For this definition of incipience, the thermal boundary -layer thickness will 
just occasionally, on a statistical basis, reach the required thickness for bubble growth 
and detachment. 

Incipience of nucleation in saturated water . - The data, listed in table I, for the 
surface- to bulk-temperature difference (AT)j at which steady emission occurred 
according to the first definition of incipience are plotted as a function of site radius 
in figure 3(a). The (AT)^ from the second definition are given in figure 3(b). The 
points presented in figure 3(a) are those (AT)^ values at which emission appeared to be 
steady but for which the next step decrease of about 10 percent in heat flux caused 
extinction or an extreme decrease of the rate of bubble production from the site. In 
figure 3(b), the (AT)^ were taken when just a very few bubbles were being produced. 
Where a range is shown about a data point, the lower limit produced no detaching bubbles 
while the upper limit provided a few bursts of bubbles on the order of 10 to 20 seconds 
apart. The incipience temperatures obtained when increasing heat flux generally were 
higher in value for a given site radius and had more scatter. The values obtained during 
increasing heat flux are not shown. The fluid bulk was maintained at a saturated boiling 
condition for all the data in figure 3. 

For a specific individual site, a range in the emission behavior would be expected as 
is characteristic of boiling processes. As shown in figure 3, the range was of the order 
of several degrees. At each site radius, the minimum (AT)j^ may be regarded as 
corresponding to conditions most favorable to bubble growth. When making comparisons 
with theory, the envelope through the minimum (AT), points should be kept in mind 
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since it delineates the boundary between the inactive region and the region of possible 
activity. 

Comparison of saturated data with analyses of incipience . - The analyses of Han 
and Griffith (ref. 4), Hsu (ref. 3), and Griffith and Wallis (ref. 5) predict respectively 
the following relations for the range of active cavity radius as a fvinction of temper- 
ature of the heated surface T . From reference 4, 



From reference 3, 



where the constants and Cg are functions of the bubble -surface contact angle and 
the angle between the surface and the tangent to the cavity mouth. From reference 5, 



The most simple of these as given in equation (4) and a modified form as given in 
reference 6 have the possible shortcoming of predicting no maximum active cavity size. 
A plot of equation (4) in figures 3(a) and (b) shows that it does lie to the left of all the 
data, but the (AT)j values are smaller than those observed. At large site radii, the 
deviation becomes larger since the data tend toward larger (AT), while equation (4) 
predicts a decrease of (AT), with increasing site radius. 

The analyses of Hsu (ref. 3) and of Han and Griffith (ref. 4) require a knowledge of 
the thermal-layer thickness. For saturated water, this thickness has been experimen- 
tally measured by Marcus and Dropkin (ref. 15) who provide a correlation between 
thermal-layer thickness and heat-transfer coefficient. For values of heat -transfer 
coefficient in the range of the present experiments (350 to 650 Btu/(hr)(ft^)(°F)), an 
average thermal-layer thickness is 6 « 0. 017 inch. 

The relation of Han and Griffith (ref. 4) as given in equation (2) is plotted in 
figure 3(a) using the thermal-layer thickness of 0. 017 inch, and it appears to bound most 
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of the data. A detailed examination of the agreement with the data, however, reveals 
some important discrepancies. For the saturated case, the maximum value of equa- 
tion (2), which is found as (AT). - ■», is R„ „„„ = 26/3 = 0. Oil inch. The theory does 
not account for the observed activity of the largest drilled sites used here which had 
R = 0. 032 inch. Possibly even larger sites would produce bubbles, so the data do not 
confirm the theoretical prediction that there is a maximum cavity radius that is active. 
Although the theoretical curve does encompass much of the data, it does not really pass 
through the lower bounds of the experimental (AT)^ points; rather it predicts (AT)^ 
values that are generally much too low. This type of comparison between data and 
theory is also evident for the measurements reported in reference 11. 

In addition to requiring a value of 6, equation (3) of Hsu (ref. 3) has coefficients 
Cj and Cg that depend on the cavity geometry. If the values of Cj and Cg used in 
reference 3 are assumed here, a curve similar to that from equation (2) is obtained. If 
as an extreme case the mouth of the cavity is assumed to have a perfectly square edge 
which would be approached by the machine drilled sites, then from the results of refer- 
ence 3 it can be shown that Cj and Cg are equal to one. With these values, equa- 
tion (3) was evaluated and is shown in figures 3(a) and (b). The curve, which strictly 
should be compared only with the mechanically drilled data because of the choice of 
constants, encompasses a little higher range of cavity radii than the Han and Griffith 
relation, but it still does not bound the data for the largest sites studied here. 

It could be argued that the nucleation that appears to come from the large cavities 
is occurring at smaller nucleation sites contained within the larger sites. However, 
Bankoff (ref. 12) shows that rectangular grooves up to a width of 0.07 centimeter 
(0. 028 in. ) will retain any vapor originally trapped within them when immersed. He 
further states that a cylindrical cavity of this (R^ = 0. 014 in. ) or somewhat larger 
radius should also retain entrapped vapor. The high-speed motion pictures, discussed 
in succeeding sections, showed a vapor cap left over the site by each departing bubble. 
Hence, it appears that the prepared sites remained filled with vapor, which makes the 
’’site within a site" explanation doubtful. The fact that the nucleus over a large nucle- 
ation site was found to grow in a saturated liquid indicates that the maximum critical 
cavity radius predicted by the analyses of references 3 and 4 may be unrealistic at 
least for the saturated case. 

The activation theories were founded on the assumption that for a bubble to begin to 
grow, the temperature in the thermal layer at a characteristic position, for example, 
at the top of the bubble nucleus, must reach the temperature of the vapor inside the 
nucleus. More realistically the growth would seem to depend on whether there is a net 
vapor addition to the bubble which in turn depends on a net heat addition to the bubble. 

If a bubble nucleus extends out beyond the thermal layer, the theory in references 3 
and 4 would predict the site to be inactive. However, the present authors believe that 
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if there is sufficient heat addition to the bubble nucleus in the region near its base, the 
nucleus may still grow even though it extends beyond the thermal layer. 

The present authors have devised a cavity activation criterion based on a heat 
balance to the bubble nucleus, and the analysis is given in appendix B. Although forming 
a heat balance appears to be a reasonable approach, the predictions were not much more 
successful than the previous analyses mentioned. As shown in figure 3(b), the results 
were close to those of equation (4) and provided values of (AT)^ that were too small. The 
analysis is included in this report for possible interest to those who may perform further 
investigations in this area. 

A summarization is as follows: Values of (AT)^ to initiate bubble production have 
been determined experimentally for a range of cavity radii and have been compared 
with theoretical predictions of cavity activation. The presently available analyses do 
not seem satisfactory for predicting the experimental observations. The observed (AT)^ 
were always significantly higher than those predicted. 


Characteristics of Bubbles at Drilled Sites 

Growth of bubbles while attached to surface . - If mechanically drilled sites in a 
polished surface and low heat fluxes near incipience were used, it was possible to grow 
isolated bubbles which were not disturbed by any adjacent bubble formation. The high- 
speed motion pictures showed the bubble profiles to be symmetric about the site axis, 
and the bubbles had little distortion. Figures 4(a) and (b) show a set of tracings for the 
growth of two bubbles, one for each of the mechanically drilled site diameters 0. 0075 
and 0. 0413 inch. These sites along with two other sites were drilled into the same 
heated strip, and the photographs were taken at the four sites in rapid succession without 
changing the experimental conditions; thus, the heat flux and the surface temperature 
were held constant. Hence, the pictures presumably compare the growth of bubbles 
under the same condition of the thermal layer on the surface. The random convecti'^^e 
circulation currents in the boiler, however, caused, some fluctuations in the thermal 
layer thickness. These fluctuations resulted in some of the bubbles at a given site having 
lower growth rates than others even though the time average conditions remained fixed. 
This statistical variation was also noted in reference 9. The statistical randomness 
must be removed as much as possible in order to help provide a consistent comparison 
of bubbles taken at different sites. The minimization of randomness was done by 
selecting the bubble at each site that had the most rapid growth rate and increased in 
volume continuously to detach in the shortest growth period. The conditions for the 
growth of such bubbles should correspond most closely to the conditions for which the 
theoretical growth equations apply. These are the bubbles shown for two of the sites in 
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figure 4 and are presumed to be growing imder thermal conditions as consistent with 
each other as possible. The surface- to bulk-temperature difference was 12. 5® F. 

The bubbles were assumed to be rotationally symmetric about each site axis, and 
their volumes were computed from the profiles. For some of the bubbles near detach- 
ment, where they were not close to spherical, the width of the bubble was measured at 
15 positions along the bubble height, and the volume was computed by numerical integra- 
tion. For the remainder of the bubbles, the shape could be approximated within an 
error of a few percent by a sphere or a spherical segment, and the volume computed by 
using spherical volume formulae. 

The volumes are given as a fvmction of time in figure 5 for the entire growth period. 
Since the photographs were taken at a film speed of approximately 2000 frames per 
second, there is a time lapse of somewhat less than 1/2000 second between frames as 
recording the image on the film requires a portion of the time. Hence, the first datum 
point can be, at most, a few ten thousandths of a second after the detachment of the 
previous bubble. If the data could be extrapolated to zero time, the initial bubble volume 
would not be zero because there was always a vapor nucleus left behind by the previous 
bubble. The last few pictures in each set of bubble profiles in figure 4 illustrate how 
the bubble neck narrows and then breaks off thereby leaving a vapor nucleus behind. 

The present data are compared in figure 5 with two of the commonly quoted simple 
equations for bubble growth: the Fritz-Ende expression (ref. 16) 



and the Plessett-Zwick expression (ref. 17) 



These expressions do not fit the slope of the data over the entire growth period, as the 
observed bubbles do not have a V variation throughout growth. The Fritz-Ende 

result, which is based on a less refined conduction model than the Plessett-Zwick equa- 
tion, agrees more closely with the range of data although the theory predicts a more 
rapid growth than that observed. The better agreement of the Fritz-Ende relation with 
the measured growth data was also found for the data in reference 18. 

In figure 6(a), the growth curves of bubbles from an artificial site of diameter 
= 0. 0075 inch are shown as a fimction of the surface- to bulk-temperature difference. 
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The curves shown were selected, as in figure 5, as being the ones which had the most 
rapid bubble growth. It was felt that these curves should correspond most closely with the 
the analytical predictions, as they presumably are for very favorable growth conditions. 
Comparison is made with the Fritz -Ende equation, and, although the experimental data 
do not fall on the analytical lines, agreement insofar as trends with AT and slopes with 
time is good. As shown by equations (5) and (6), the bubble volume should depend on 

O 

(AT) , and this relation is indicated to some extent by comparing the data of figure 6(a) 
to the curves of the Fritz -Ende relation. The variation of volume with AT will be 
discussed in more detail in connection with figure 7. 

In figure 6(b), growth curves for bubbles originating at a larger site with diameter 
= 0. 0413 inch are shown. The correlation with the analytical curves, and the self- 
consistency of the data is obviously poorer than for the smaller site. Two possible 
reasons for poorer correlation are evident. First, the models used in deriving the 
analytical expressions have the assumption that the bubble is surrounded by superheated 
liquid. This assumption is of doubtful validity for the very large bubbles growing from 
the larger sites especially during the later stages of growth. Also, because of the lower 
frequency of bubble emission from the larger sites at a given AT, relatively fewer 
bubbles were available on each film roll for obtaining growth curves. Thus, in selecting 
the most rapid growth curve for each given AT, the probability of obtaining comparable 
curves for comparison between the different AT values was not as good as for the higher 
frequency small site. 

In figure 7, the growth data of figure 6(a) for bubbles nucleated from the 0. 0075- 
inch diameter site are cross -plotted to show more clearly the relation between bubble 
volume and AT at fixed times from the beginning of bubble growth. The increase in 
bubble volume with (AT) as predicted by the Fritz-Ende equation is shown by the dotted 
lines, and agreement with the limited data available seems reasonable. 

Interpretation of bubble-growth results . - Some aspects of bubble departure from a 
known site have been studied previously by Wei and Preckshot (ref. 19), who photo- 
graphed bubbles leaving a short section of glass capillary that had been cemented to a 
glass slide placed on a heated copper block. The entire assembly was submerged in 
water and bubbles nucleated from the capillary. In this case, where the cavity itself 
was not within the heat generating material, after a bubble detached there was a small 
vapor pocket left at the bottom of the cavity, and the remainder of the cavity filled with 
liquid. The small vapor pocket grew to form the next bubble. In the present experi- 
ments, the cavity was drilled into the heat generating material of the test strip, and 
either because of this or the cavity shape or material, the vapor nucleus was much 
larger than that in reference 19. The vapor nucleus was always observed to protrude 
well above the cavity opening. For this reason, it is possible that the walls within the 
cavity were dry without any liquid layer penetrating inside. Within the very short time 
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after bubble departure and before the next frame showing the new bubble, the surface 
tension formed the nucleus into the shape of a spherical surface. For a small site 
(fig. 4(a)), the nucleus was almost in the shape of a complete sphere, while for a larger 
site (fig. 4(b)), it was approximately a spherical segment. 

It is postulated herein as a possible growth mechanism that the growth of the bubble 
depends on the heat addition from the thermal layer at the bubble base. Hence, the 
growth of the bubble depends on how much surface area it has within the thermal layer. 
Analyses such as those by Zuber (ref. 20) and Han and Griffith (ref. 4) assume that the 
growing bubble nucleus pushes the thermal layer outward; thus the bubble is surrounded 
by the thermal layer. As was discussed in the section Growth of Bubbles While Attached 
to Surface, the thermal-layer thickness for the heat fluxes encountered here was on the 
order of 0. 017 inch in thickness. This thickness is shown in figure 4 and is a little 
larger than the height of the vapor nucleus left behind by the departing bubble. During 
most of the bubble-growth period, the thermal layer was small compared with the 
bubble height. This comparison was true for the cases reported herein because the 
artificial sites were fairly large, and large bubbles were produced. As a result, it is 
possible that the heat transfer to a bubble was only through the area of the base in 
contact with the thermal layer. This area is identified as the area A in figure 8. The 
growth of the bubble then depends on the time variation of A. As a very simple illustra- 
tion, assume that the heat -transfer coefficient at the bubble base and the average tem- 
perature within the thermal layer remain constant throughout bubble growth. Then, the 
energy added from the liquid to the bubble is equated to the energy required to provide 
vapor for bubble growth: 


hA(Tj^ -T^) = ^p^X (7) 

dr 

The time variation of A then must be considered in order to integrate for V. From the 
bubble profiles in figure 4, after about the first 0. 001 second the bubbles are quite close 
to spherical segments and are of the successive configurations shown in figure 8(b). The 
area A of any spherical zone is given by ttDL, where L is the height of the zone. 
Therefore, the area of the spherical zone within the thermal layer is A = t 7D6. If the 
volume of the bubble is approximated by V = ;rD^/6, then A = n-(6V/jr)^^^6. Substituting 
into equation (7) and assuming that the thermal-layer temperature and thickness remain 
constant during bubble growth give the functional form 

yl/3 „ W 

dr 
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or 


V oc (8) 

In the later stages of growth when the bubble begins to elongate and especially when 
a neck begins to form, the area A tends to become constant, especially for the bubbles 
from large sites as in figure 8(c). Then from equation (7), dV/dr = constant, which 
gives 


V oc T (9) 

This first power behavior of the volume with time was also observed in reference 18. 

The trend of decreasing slope during growth of the curves in figures 5 and 6 then might 
be attributed in part to the changing configuration near the bubble base. The fact that 
the temperature in the thermal layer T^^ actually decreases during growth (refs. 21 
and 22) and that the heat added from the thermal-layer thickness 6 may also decrease 
during growth (ref. 23) contributes further to a reduction in growth rate. 

When a bubble is formed at a very small natural nucleation site, the edge of the 
bubble base grows beyond the limits of the site opening. The base expands rapidly out- 
ward at first, and then the base diameter remains relatively constant imtil bubble 
detachment begins (e. g. , base-diameter measurements in ref. 18). The rapid growth 
in the base diameter produces a rapid increase in the bubble area within the thermal 
layer which could accoimt for the fact that in the expression (ref. 24) 

D cc 7^ (10) 

the n values have been as large as one near the beginning of bubble growth. 

Forces acting on bubbles at departure . - The shapes of typical bubbles formed at 
drilled sites are shown in figure 9 just prior to the beginning of necking at the base of 
the bubble, which is the start of bubble departure. In the present experiments, the 
forces on the bubble at departure are believed to be computed very accurately because 
of the following considerations: At departure, the cylindrical neck at the base of the 
bubble was found to stretch in length and become narrower imtil the neck was broken. 
Hence, the force holding the bubble to the surface is the surface tension of the interface 
around the circumference of the bubble neck. Consequently, for the bubbles considered 
herein, the surface-tension force can be computed quite accurately since it does not 
depend on the contact angle where the liquid-vapor interface meets the heated surface. 
During boiling, this contact angle is difficult to measure experimentally with high 
accuracy. As will be shown, the diameter of the bubble neck is known accurately here 
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because it corresponds quite well with the diameter of the drilled site. 

The buoyancy force of the bubble at departure can be computed by a method proposed 
in reference 25. The bubble parameters that were obtained from the photographs are the 
total volume V, horizontal maximum diameter D, height H, and base diameter D|^ 
as shown in figure 9(d). The buoyancy expression will be placed in terms of these 
quantities. The buoyancy of the bubble volume, excluding the shaded volume V^j of 
figure 9(d), is given by 


(Pi - Py) -S- (V - Vl,) 
®0 


Assume the shaded volume is approximated closely by a cylinder (thus, neglect the 
curvature of the top of the bubble); this term can then be written as 



( 11 ) 


The contribution of to the buoyancy force can be expressed in terms of the pressure 
difference across the interface at the top of the bubble 



(Py - Po> 


( 12 ) 


This pressure difference can be written in terms of the two principal radii of curvature 
at the top of the bubble 


(Pv - Po> = ^ 




top of bubble 


(13) 


For the bubbles given herein, the upper portions are so spherical that and R£ at 
the top of the bubble can both be approximated by Rj^ = Rg = D/2. The sum of equa- 
tions (11) and (12) with the use of equation (13) provides the final expression used to 
compute the buoyant departure force. 


Py> 


sA 4 / D 


(14) 
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As discussed earlier, the bubble departs by breaking the neck a little above the 
position where it is in contact with the surface. The surface-tension force is then given 
by 


F,p = ottDjj (15) 

The width of the neck at the bubble base was measured from the photographs. As shown 
by the values in figure 9, the width was found to be quite close to the diameter of the 
drilled site except for the smallest site photographed where in one instance the neck 
appeared somewhat wider than the site. Since the drilled holes were measured directly 
with a microscope, the precision is higher than the photographic measurement at the 
bubble base as the photographs are distorted somewhat by density gradients in the fluid. 
Hence, using the diameter of the drilled site and computing the surface tension force 
from 


= (nrD^ (16) 

is considered more accurate. As shown in the previous section on bubble growth rates, 
the volume change with time for the present bubbles was small at the time of departure. 
This small growth rate eliminates the possibility of inertial forces being of significance. 
It also eliminates the consideration of drag forces resulting from the motion of the ex- 
panding surface of the bubble. 

A summary of bubble -departure dimensions and forces is given in table n. A plot 
of buoyancy against surface-tension force is given in figure 10. Within a statistical 
scatter that seems reasonable for a boiling experiment, the agreement of the two forces 
is very good. This agreement shows that for the large slowly growing bubbles studied 
herein, the buoyancy force defined as the integrated vertical pressure force on the bubble 
is the force causing departure. 

Equilibrium bubble shapes . - In reference 26, Bashforth and Adams studied the 
shapes of drops and bubbles at equilibrium subject to the forces of buoyancy and surface 
tension. They provide tables of coordinates of the bubble profiles in terms of the 
parameter 


g(Pv - P;)b^ 

/3 = — 1 i 

go<^ 

Since some of the present bubbles fall so close to the correlating line in figure 10, which 
indicates an equilibrium between surface tension and buoyancy forces at departure, it is 
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of interest to see whether these bubbles compare with the predicted equilibrium shapes. 
Only the shapes are being compared and not the actual bubble sizes since the theoretical 
results are in terms of dimensionless heights and radii utilizing the radius of curvature b 
at the top of the bubble as a scale factor. Multiplying by the dimension b to place the 
curves in dimensional form would not change the relative shapes at a fixed /3. Since 
is greater than p^, the values of ^ are negative for the present results. The two 
typical experimental bubbles shown in figure 11 have values of 0 = -0. 24 and -0. 35 and 
grew at mechanically drilled site diameters of 0. 0209 and 0. 0413 inch, respectively. 

The tables of bubble coordinates in reference 26 are somewhat incomplete in the range of 
/3 = -0.1 to -0. 3, and the theoretical curves were supplemented with values given in 
reference 27. 

The comparisons in figure 11 reveal that the upper portions of the experimental 
bubbles are very close to spherical and agree well with the theory. At the base portion 
of the bubbles, the experimental profiles deviate from theory probably because of the 
physical constraint imposed by the attachment of the bubble base to the opening of the 
drilled site. The theory on the other hand does not have this constraint. These con- 
siderations reveal an additional factor, other than dynamic forces present for rapidly 
growing bubbles, which would account for why the Fritz equation (ref. 28) does not always 
correctly predict the size of bubbles at departure. The Fritz equation is based on theory 
which assumes that the base of the bubble is not constrained while, depending on the site 
diameter and geometry, the base of some actual boiling bubbles can remain attached to 
the site as in the present study. 


CONCLUSIONS 

1. Bubble incipience is a vague concept in relation to interpretation of experimental 
data. Analytical determinations of site activation consider the surface superheat 
required to make a bubble nucleus grow although the bubble may not then detach. Incip- 
ience is defined here in terms of the condition when a steady bubble stream is being 
formed. This condition requires that the bubble must not only begin to grow, but it 
must also be able to reach a size sufficient to cause detachment. 

2. Experimental values were obtained for the surface superheat when a cavity begins 
to produce bubbles, and the results were compared with several analyses. For the 
range of cavity radii tested, 0.0015 to 0.032 inch, the predicted (AT)^ were always 
significantly smaller than those measured. 

3. The experimental results show that there can be active cavities which are larger 
than those indicated by some of the analytical models relating cavity size and thermal- 
layer thickness for site activation. These models predict that the maximum active site 
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radius is comparable to the thermal-layer thickness, whereas experimentally much 
larger sites were found active. This conclusion is based on the assumption that a ther- 
mal layer of thickness 0. 017 inch is present in the present experiments. The upper 
limit of active site radius is probably governed in part by the ability of liquid to fill the 
site rather than by considerations of thermal effects. 

4. In the initial period of bubble growth, the bubble volume was found proportional 
to (time)^/^. Later in the growth period, the volume is proportional to time, which 
implies that vapor addition to the bubble may take place through a constant area around 
the bubble neck. 

5. For a given small site, the volumetric growth rate increases approximately as 

O 

(at) as indicated by many of the growth equations, where AT is the difference between 
the temperature at the boiling surface and the saturation temperature. For a large site, 
the functional relation is not clear, but the growth rate does increase with AT. 

6. The bubbles detach from the artificial sites studied here by forming a neck which 
is ruptured some distance above the site. Consequently, a vapor nucleus is left behind 
from which the next bubble grows. 

7. For slowly growing bubbles, the buoyancy and surface-tension forces are in 
balance at bubble detachment. The buoyancy force must include the pressure force def- 
icit caused by the missing bubble surface area in the contact region at the bubble base. 

8. Although uniform temperature conditions were not present, the bubbles near their 
departure from the artificial sites compared reasonably well with the equilibrium shapes 
predicted by Bashforth and Adams. Near the bubble base, the constraint imposed by 
attachment to the artificial site causes some deviation from the predicted shapes. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 13, 1967, 

129-01-11-06-22. 



APPENDIX A 


SYMBOLS 


A 

surface area 

I^c 

radius of artificial nucleation 

a 

thickness of specimen contain- 


cavity 


ing artificial site 

Rl , R 2 

principal radii of curvature of 

b 

radius of curvature at top of 


bubble 


bubble 

AT 

temperature difference. 

Cl, Cg 

dimensionless site geometry 


T - T 
w sat 


constants 

"^bulk 

bulk temperature of boiling 

D 

maximum horizontal bubble 


liquid 


diameter 

To 

temperature at insulated sur- 


base diameter of bubble 


face of test specimen 


diameter of artificial nucleation 
cavity 

Tsat 

saturation temperature of 
boiling liquid 

Fb 

Ft 

buoyancy force 
surface-tension force 

TtZ 

Tv 

mean thermal-layer temper- 
ature 

temperature of vapor inside 

g 

gravitational acceleration 
conversion constant, 32.2 

V 

bubble 

®0 

H 

h 

Ib-mass-ft/lb force-sec^ 
bubble height 

boiling heat -transfer coefficient 

T 

w 

V 

Vb 

temperature of boiling surface 
of test specimen 

volume of bubble 

bubble volume which lies above 


evaporating and condensing heat- 
transfer coefficients for heat 
transfer to bubble nucleus 

X 

area of bubble base 

distance coordinate from sur- 
face 

k,k^ 

thermal conductivity of test 
specimen or liquid 

a 

thermal diffusivity 
bubble -shape parameter. 

Pv’Po 

internal and external pressure 
at top of bubble 


g (P, - P;)b2 

Q 

energy per unit time 


So 

q 

heat flux 

6 

thermal-layer thickness 

R 

20 

bubble radius 

X 

heat of vaporization 


Subscripts: 


density of liquid or vapor 

surface tension of liquid-vapor 
interface 

time 


bulk 

i 

sat 

V 


fluid bulk condition 
condition at incipience 
saturated 


vapor 


APPENDIX B 


CAVITY ACTIVATION CRITERION 

The analysis presented herein is based on an energy balance for a bubble nucleus that 
is assumed to exist at a nucleation site. As discussed in reference 5, as a vapor nucleus 
grows out of a surface cavity, the minimum radius of curvature of the vapor -liquid 
interface is reached when the vapor cap becomes a hemisphere. The hemispherical 
nucleus will be the most difficult bubble configuration to grow and will thus have the 
hipest vapor temperature within it. For this reason, the analysis will deal with the 
growth of a hemispherical nucleus. This nucleus is represented in figure 12. Although 
the data in the present report are for saturated bulk liquid, the analysis will include the 
subcooled case. 


Temperature Profile 


As shown in figure 13, for simplicity the temperature distribution in the liquid 
thermal layer near the surface will be approximated by a strai^t line extending from the 
wall to the thickness where the bulk temperature is reached. The thickness 6 

designates the position at which the fluid is at the saturation temperature. The temper- 
ature profile is then given as 


■ '^bulk ^"^w 



X ^ 6 


bulk 


(Bl) 


TW = T^ulk 


X 5: 6 


bulk 


(B2) 


Growth Criteria for Bubble Nucleus 


A hemispherical vapor nucleus is now assumed to be on the surface cavity. This 
assumption is somewhat restrictive since for high subcooling, especially in conjimction 
with a large cavity size, the nucleus may not be able to grow this large. The vapor 
temperature within the nucleus is assumed to be constant and is given by (ref. 5) 


2oT 


Ty - Tsat 


sat 


Ap R 
^v c 


(B3) 
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Three cases can be distinquished depending on the relation of cavity radius to the 
thermal -layer thickness. One case is shown in figure 12(a) where the nucleus extends 
through the thermal layer (R^ 2 : Figure 12(b) where < ^bulk the 

other two cases depending on the size of R^ relative to 6^, the position at which 
T(x) = T^. The two cases are thus for ^ 6^ and for 0 ^ R^, < 6^. If equa- 

tion (B3) is used for the vapor temperature within the bubble, the net heat addition to 
the nucleus is 




+ h^27TRj. 




"^bulk ^^w " '^bulkH^ 


^bulk/ 


V 

2orT 


- Tsat - l<^ «bulk (B«) 




Q = \2,R^ 


ft 


Tu.,n, + (T - Tv 


bulk ' w ” bulk^ 


) l-_iL_ -T 


2oT 


‘'bulk 


sat 


sat 


^v^c 


dx 


0 ^ Rc - ^ 


(B4c) 
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The quantities h and h are the evaporating and condensing heat-transfer coefficients 
for heat transfer to the nucleus and have been assumed constant. For the nucleus to 
grow, Q must be greater than zero. Solving equations (B4) for the activation criterion 
Q > 0 and assuming h. = h result in the following two relations 

(^w ■ ^sat). ~ ~ ^"^sat " ^bulk^^^ ^c “ ^bulk 

^ ^v^bulk '^®bulk > 


(AT)i 


2oT 


(T 


sat ” "^bulk^ 


sat 


25 


bulk 


APyRc 1 - 


R, 


25 


1 - 


R, 


^c “ ®bulk 


(B5b) 


bulky 


25 


bulk 


Equation (B5b) incorporates the incipience criterion for both (B4b) and (B4c). For 
saturated conditions “ "^gat ®bulk ~ (B5a) and (B5b) become 


sat ^ 


4oT 


sat 




Rc^6 


(B6a) 


sat 



R_ :£ 6 


(B6b) 


\Vhen equations (B5a) and (B5b) are compared with equations (B6a) and (B6b) it is evident 
that relative to the saturated condition the subcooling provides an additional term that 
increases the temperature difference required for growth of the nucleus. No maximum 
site radius that will be active is indicated by this analysis. 

A similar analysis except that it is based on the surface heat flux rather than the 
surface temperature has been given by Grace (ref. 29). 


Temperature Difference Required for Bubble Detachment 

For a saturated bulk fluid, if the activation AT .is exceeded, the theory predicts 
that bubbles will grow until they detach from the surface. When the bulk is subcooled, 
however, the upper portion of the bubble may protrude into liquid that is below the vapor 
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temperature within the bubble, and condensation will occur. As a bubble grows, a size 
may be reached where the condensation balances the evaporation in the region near the 
bubble base. Then the bubble can no longer increase in size. If this zero net heat 
addition condition occurs before the bubble has reached the size necessary for detach- 
ment, the bubble will not leave the surface. Then the site is not incipient because no 
stream of bubbles is being produced. Thus, in the case of a subcooled fluid, two 
criteria have to be examined with regard to the definition of incipience. The first is 
whether a vapor nucleus at the nucleation site will grow; the second is whether a bubble 
can become large enough to detach from the surface and produce a stream of bubbles. 

To arrive at a criterion for whether bubbles will grow to detachment, the size a bub- 
ble will achieve must be expressed in terms of the subcooling. For simplicity, the bub- 
ble is assumed to be a portion of a sphere as shown in figure 14. The bubble radius, 
height, and the cavity radius are then related by 


R = 





2H 


(B7) 


The final expressions may be simplified by the assumption that the evaporation and the 
condensation heat -transfer coefficients are eqiial. Then the heat addition to the bubble is 


Q = 27rRh 


If 


'bulk 




('^w ■ "^bulk^ " 7^ I * '^bulk " '^sat 


2oT 


sat 


XP^R 


dx 


Q = 27rRh 


/'(’ 

^bulk 


2oT 


'^bulk ■ "^sat 


sat 




|dx > H > (B8a) 


f 




^bulk^ 


2oT 


*^bulk ■ "^sat 


sat 


XP^R 


* 'bulk 


The heat addition Q is set equal to zero to provide the condition when the bubble will 
remain of fixed size. The integrals are carried out, and the radius R is eliminated in 
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terms of the height H by using equation (B7). The result is an implicit expression for 
the height of the equilibrium bubble, 


H ^ i/ '^w " ^sat 
^buUc ^ V sat " "^bulk 


+ 1 - 


4oT 


sat 




H > 6 


^v^bulk^'^sat " "^bulk^ + R? 


bulk 


(B9a) 


H 

^bulk 


2(Tw - Tsat) 


<^w - Tsat) + (T 


sat ” ^bulk^ 


8orT 


sat 


H 




(T^ - Tgat^ ^'^sat " ^bulk^] 


H < 6 


bulk 


(B9b) 


If a given set of variables such as subcooling, cavity radius, etc. , is inserted into equa- 
tion (B9), the equilibrium bubble height can be found, and the question then remains as 
to whether a bubble of this height will depart from the surface. 

The departure depends on the balance of forces on the bubble. The most simple 
case, which was valid for the large sites in saturated liquid studied herein, is when the 
departure depends only on the balance of buoyancy and surface tension. By equating 
equations (14) and (16), the bubble height at departure is obtained in terms of the cavity 
radius. The height can then be compared with that from equation (B9) to determine 
whether departure will occur. 
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TABLE I. - TEMPERATURE DIFFERENCE BETWEEN SURFACE AND SATURATED 


LIQUID FOR STEADY BUBBLE PRODUCTION, SATURATED CONDITIONS 


Electron -beam 

Machine drilled sites 

Electron -beam 

Machine drilled sites 

drilled sites 



drilled sites 



Site 

Surface- to 

Site 

Surface- to 

Site 

Surface- to 

Site 

Surface- to 

radius, 
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TABLE n. - BUBBLE DIMENSIONS AND FORCES AT ONSET 


OF BUBBLE DEPARTURE IN SATURATED WATER 


[Mechanically drilled sites.] 
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Cratered lip 


(b) Obtained by mechanical drill; depth, 0.010 
Inch. 290X. 


(a) Obtained by electron-beam disintegration 
depth, 0.0025 inch. 250X. 


(c) Obtained by mechanical drill; depth, 0.020 Inch. 150X, 


re 1. - Artificial nucleatlon sites. (All dimensions In inches.) 




Insulated stain less - 
steel tank 



Figure 2. - Schematic diagram of experimental apparatus. 
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Equation 

Thermal-layer 

thickness, 

6 , 


in. 

(4) 

(2) 

(3) (Cl = C 3 ' 1) 
(B 6 ) 

0.017 

.017 

.017 

Drilling method 


Electron beam 
Machine 
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□ 
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T\ “0TT" 


(a) Incipience temperature based on regular site activity. 






hCH KIH 


I 2 4 6 8 10 12 14 1( 

Surface- to saturation -temperature difference, (AT)j, °F 

(b) Incipience temperature range from where there are few bursts of bubbles to where bub- 
ble detachment ceases. 

Figure 3. - Temperature difference required to initiate bubbles from artificial sites in water 
at saturation temperature. 







Time, T, sec 0.000455 . 00137 . 00318 . 00591 

Volume, V, in. ^ 5.65x10*^ 4.00x10'^ 1.47x10'^ 3.85xl0’5 



Time, t, sec 
Volume, V, in. ^ 


Time, t, sec 
Volume, V, in.^ 



0.0300 

2.64x10"'’ 


(a) Site diameter, 0.0075 inch. 



Volume, V, in.-^ 


9. 18x10 9. 30x10“^ 

(b) Site diameter, 0.0413 inch. 


Figure 4. - Growth profiles of bubbles at artificial site. Temperature 
difference AT, 12. 5° F; thermal layer thickness 6, 0.017 inch; 
heat flux q, 7100 Btu per hour per square foot. 





Bubbie volume, 








Growth time, 



ence, AT, °F 

Figure 7. - Effect of surface- to saturation- 
temperature difference on bubble growth 
rates from artificial sites. Site diameter, 
0.0075 inch. 
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site diameter, 0.0209 inch site diameter, 0.0413 inch 

(from microscope). itrom microscope). 

Figure 9. - Typical bubble configurations just prior to necking and departure. 
Temperature difference, AT, 12. 5° F. (All dimensions in inches. ) 
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igure 11. - Comparison of theoretical and experi- 
mental bubble profiles at onset of bubble departure 
when buoyancy and surface tension forces are in 
equilibrium. Temperature difference, AT, 

12.5° F. 



Temperature 

Figure 13. - Significant temperatures and thicknesses in thermal 
layer. 



Figure 14. - Truncated sphere 
bubble model. 
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*'The aeronautical and space activities of the United States shall be 
conducted so as to contribute to the expansion of human knowU 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof/* 

— ^National Aeronautics and Space Act of 1958 
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